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Abstract 

We consider a size-structured model for cell division and address the question of determining the 
division (birth) rate from the measured stable size distribution of the population. We formulate 
such question as an inverse problem for an integro-differential equation posed on the half line. 
We develop firstly a regular dependency theory for the solution in terms of the coefficients and, 
secondly, a novel regularization technique for tackling this inverse problem which takes into account 
the specific nature of the equation. Our results rely also on generalized relative entropy estimates 
and related Poincare inequalities. 

Key-words Size structure models. Inverse problem. Regularization technique. Entropy methods. 
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1 Introduction 

For many unicellular organisms the mass of the cell, its DNA content, and the level of certain proteins 
concentrations are often considered as the most relevant parameters for modeling the cell division and, 
consequently, the population dynamics. Thus, it is natural and usual to consider the evolution of the 
cell density n as a function of the time t and a size parameter x. See [H EJ [3] for applications to 
biology and [I] for other related cases. In many situations it is of paramount importance to determine 
from observed data the division (birth) rate B as a function of the size parameter x. However, one is 
faced with the problem that the cell distribution n(t, x) is not easily measured or observed. Instead 
of that, the dynamics of n(t, x) leads to a stable steady distribution N = N(x), which after a suitable 
re-normalization, can be observed. In this work we face the problem of determining B = B{x) from 
noisy and scarce observed data N = N{x). 

Age-structured populations also constitute basic models for the division of certain cells, [H[6]- There 
has been substantial interest on the inverse problem for structured population models. See for example 
EJ El Uni E] and references therein. These models typically lead to difficult inverse problems that 
require a substantial amount of regularization, especially if one is interested in making use of real data. 
In full generality, size structured models have however a somewhat different mathematical structure 
(only under the assumption the smallest daughter cell is larger than half the largest mother cell, the 
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models lead to similar theory). For this reason, our approach, differently from that of other works we 
are aware of, makes use of the information given by the so-called stable size distribution of the model 
and relies on a number of recent methods based on generalized relative entropy estimates. 
As a starting point, we consider the following size-structured model: 



Jjn(t, x) + J^[g(x)n(t, x)\ + B(x)n(t, x) = AB(2x)n(t, 2x), x > 0, t > 0, 
n(t,x = 0) = 0,t > 0, 
n(0,x) = n°(x) > 0. 



(1) 



This expresses that the evolution of the cell density results from two effects: At the one hand, the 



term Q^[g(x)n(t,x)] which describes the growth of cells by nutrient uptake with the rate g(x) 



At 



the other hand, the terms containing B(x) and B(2x), which describe the division of cells of size 2x 
into two cells of size x. The resulting function-differential equation has a mathematical interest per 
se. One way to understand such interest is that the x-derivative leads to transport to right, whereas 
the functional term leads to transport to the left. The overall effect is an equilibrium. Another way 
to understand those effects, in mathematical terms, is to consider the evolution of two macroscopic 
quantities, the total cell number N(t) = J °° n(t, x)dx and the total biomass M(t) = xn(t, x)dx. 
Integrating the equation (pQ) yields 



d 



N{t) 



B(x)n(t, x)dx, 



dt jo 

this indeed means that number of cells increase only by division, and 



d r°° 



(x)n(t, x)dx 



(2) 



(3) 



which indeed means that the biomass increases only by nutrient uptake. 

For g = 1, a case we consider here for simplicity, it was shown in [12|, I13j that there exists a unique 
eigenvalue Ao and eigenfunction N = N(x) such that, after a suitable time re-normalization, the 
solutions of dU) converge to a multiple of N thus given by the solution of the eigenvalue problem 



+ (Xo + B(x))N 



N(x = 0) = 0, 



AB(2x)N{2x), 



x > 0, 



(4) 



N(x) > for x > 0, J °° N(x)dx = 1. 
Namely, we have under fairly general conditions on the coefficients 



n(t, x) 



-Ao* 



pN(x) 



in weighted L p topologies that are related to entropy properties and that will be described later on. 
Moreover exponential rates have been proved to hold for fairly general rates B ( \\.2\ [2]). Such N is 
therefore the above mentioned stable size distribution] it is the distribution observed in practice and 
available for measurements. 



The precise question under consideration in this work is the following: 
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Inverse Problem 1.1 How to recover in a stable way the birth rate B(-) from noisy data N(-) and 
the rate Ao ? 

The plan for this work is the following: In Section [2] we present some preliminary remarks on a toy 
model that describes and motivates our approach to the regularization of the inverse problem under 
consideration. In Section [3] we collect a number of results concerning the direct problem associated to 
the model. In particular, concerning the dependence of the solutions to the model with respect to the 
coefficients. In Section |4] we address the inverse problem and its regularization. The regularization we 
propose is novel and consists of introducing the regularization parameter a directly in the equation as 
a coefficient of the perturbed differential equation. We establish a result concerning the strong stability 
of such perturbed equation and another one demonstrating the consistency of such perturbation when 
a goes to zero. In Section [5] we obtain a convergence rate for this regularized inverse problem. Some of 
the proofs of the results in Section [3] are rather technical and we post-pone them to the Appendix lAl 



2 Preliminaries: A Classical Example Revised 

In this section we present a simplified problem which shares some similarity with and allows us 
to present our strategy for attacking the Inverse Problem 11.11 This problem is the regularization of a 
classical ill-posed inverse problem, namely recovering a function from its antiderivative, from a slightly 
different perspective. We believe this will make our approach and our estimates more clear by looking 
at a specific and simpler well known example. We refer the reader to the books |14t [T5] for further 
information and the classical treatment of the subject. 

Inverse problems are characterized by the fact that they are typically ill-posed. Perhaps one of the 
best well-known examples is that of differentiation. Namely, find a function u from its antiderivative 
v. We focus on the following: 

Inverse Problem 2.1 Find u : I — > R such that 



X 



X(s)u(s)ds = v(x) , x £ I (5) 







where I Cl is an interval (finite or infinite) of the line with its left endpoint and A : I — > M>o is 
given weight function. 

Equation © is formally equivalent to 

An = °-v. (6) 

We consider the problem ([5|), and notice that the map that sends v i — ► u is not continuous, when 
we endow the space of the v's with say the L 2 (I; dx) norm and the space of the u's with L 2 (I; X(x)dx) 
norm, the choice of spaces being motivated by the specific extansion we have in mind. In practice, 
this reflects itself on the fact that small changes in the measurement of v may lead to huge changes in 
u. Still, one is forced to face such problem in a number of applications. We now describe some of the 
main ingredients from inverse problem theory that are used to handle such difficulties as they will be 
applied to the context of this article. They consist of the following ingredients: 

1. Assume some kind of "a priori" bound on the solution u to problem ([5]), or equivalently on the 
acceptable data in a sufficiently strong norm. 
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2. Assume that the measured data v e satisfies \\v — v e \ \ < e in a suitable weak norm || • ||. 

3. Approximate the problem ([5]) by one that is well-posed for the weak norm and solve the stabilized 
problem instead of the unstable one choosing the regularization parameter in an optimal way. 

Let us illustrate the technique in the problem under consideration. For more information we refer the 
reader to [HI [15] . Let us consider 

a^(Xu a ) + Xu a = £v , 

Xu a {0) = { > 

In what follows we shall assume that v(0) = whenever v is defined in a neighborhood of due to 
some extra assumptions (such as j^v € L 2 ). 

In order to implement the above program to the problem under consideration, we shall establish the 
following two inequalities: 



Lemma 2.2 If u a satisfies then 



\ u a\\L 2 {I,\ 2 dx) < -\\ v \\L 2 {I,dx) ■ ( 8 ) 



a 



If -§£iv E L 2 (I,dx) and u a satisfies th 



en 

d 2 

\Wa\\L 2 (I,\ 2 dx) < a\\-Q^2 V \\L 2 (Ldx) ■ (9) 

As a direct consequence of Equation ([9]) we have that if u a is a solution of equation ([7]) and u is a 
solution of ([6]) , then the difference u a — u satisfies 

. . d 2 M 

\\U a - -"||L 2 (/,A 2 dx) < a||^2^llL 2 (/,dx') ■ V 10 ) 

We now join the above estimates with the "a priori" assumption that 

\\^ V \\L*{I,dx)- E ' (U) 

and consider a problem with noisy data given by v e such that 

\\ v ~ v e \\ L 2(i,dx) < e • (12) 
Thus, upon solving the regularized problem © with data v e we are left with a solution u e a of 



A<(0) = . 

We further assume that J^v t G L 2 (I,dx). Now we consider the distance between u e a and u 

\\ u a - u\\L 2 (I,\ 2 dx) < \\ u a ~ u\\l 2 (I ,X 2 dx) + \\ u a ~ u a\\ L 2 (I ,X 2 dx) 



(13) 



d 2 I 
a ' 



< ®\\-^v\\ L 2(j 4x) + - \\V € - V\\ L 2(J 4X ) 



< aE+- 

a 



4 



Now, we are free to choose a in such a way as to minimize the R.H.S. of the last inequality. Using 
the elementary fact that a i— > (aE + e/a) has a minimum at a = \/E/e and choosing such optimal 
a = a(e) we get that 

\K-u\\L*a,x*dx) <2e 1 / 2 E 1 / 2 . (14) 

Summing up, if we know an "a priori" bound on the data v and v e of the form (jl II) and the data is 
measured with an accuracy e then we get that the stabilized computed solution is within 0(e 1 ^ 2 E 1 ^ 2 ) 
of the actual solution u in the norm of the space L 2 (I, \ 2 dx). 

The reasoning developed above can be pushed in different directions and in the next subsection we 
will show how it can be extended further. 

We remark that Equation (|7|) is just one of the different possibilities of regularizing the Problem 
([6]) . Another possibility would be to solve an equation of the form 

d 2 u d 

~ a ^~9 ~T + = ( 15 ) 

ox z A ox 

subject to appropriate boundary conditions. It is an easy exercise to verify that Equation (|15p formally 
corresponds to Tikhonov regularization of the Inverse Problem 12. 11 The form of Equation (|15|) . 
although very attractive and general is not well adapted to our hyperbolic problem and we do not 
consider it here. 



3 Model Properties and Regular Dependency upon the Coefficients 

In this section we shall address a preliminary issue, namely the smoothness of direct map 

B i-> (A , N) 

in suitable function spaces. Then, the Inverse Problem 11.11 makes sense in such spaces. This requires 
some assumptions and preliminary mathematical material that we introduce now. 

Firstly, we need some assumptions. Here we make the choice of simplicity, keeping in mind that 
more realistic hypotheses are possible that would lead to more technical proofs, which we prefer to 
avoid. The measurable division rates B are supposed to satisfy 



3B m , B M such that < B m < B(x) < B M , 



(16) 



Secondly, because the adjoint equation is as important as the direct eigenvalue problem, we consider 
both the cell division equation, to find (Ao, iV) solution to 



( 9 



£N + (A + B(x))N = 4B(2x)N(2x), 



x > 0, 



N(x = 0) = 0, 
N(x) > for x > 0, 



(17) 



f °°N(x)dx = l. 



and the adjoint equation 



d_ 



9?-(Ao + £(£)V = -2£(xM§), 



x > 0, 



ip(x) > 0, / °° ip(x)N(x)dx = 1. 



(18) 
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The solution <p is used in several places because it gives natural bounds for the direct problem. See 
the proof of Theorem 13.21 and the entropy dissipation in Appendix A.l. 



The existence and uniqueness of solutions to these eigenproblems was proved in [12] , and it is known 
that iV as well as its derivative vanish at x = and x = oo. We also recall the following results that 
we will use later on. 

Theorem 3.1 \12\ {Tfi^ Under the assumption Iil6\). the solution N to ( [i7| ) satisfies the properties 

B m <X = J B(x)N(x)dx < B M , (19) 

-L < I xN(x)dx = -L < (20) 

< N(x) < 2B M , (21) 
e ax N G L 1 nL°°(M+), Va < A + B m , ip(x) < C(l + x). (22) 



Proof. We only prove the estimates on N and refer to |16| for the sub-linearity result on ip and for 
more accurate results on the behavior of N at infinity, in particular, when B has a limit B^ at infinity, 
one can obtain exponential decay with a = Ao + B^. 

The inequalities (|19p and (|2U|) are obtained by integrating Equation (|17j) against respectively the 
weights dx and x dx. 

For (|21h . we use (I19h . and from the Equation (|17p . we obtain 

-0-N(x) < AB(2x)N(2x) =^ N(x) < 4 I B(2y)N(2y)dy < 2B M - 
ox Jo 

For (|22|) . we multiply (|17|) by e ax and integrate between and oo. In fact, a standard truncation 
argument is needed for complete justification. See [2] for details. We obtain 

/ °°[Ao + B(x) - a]e ax N(x)dx = 4 1 °° B(2x)N(2x)e ax dx 

(23) 

< 2/ °° B{x)N{x)e ax / 2 dx. 
With the choice a = Ao, since the right hand side is bounded by 

/ [-co \ 1/2 / .ao \ 1/2 / [-co \ 1/2 

2( / B(x)iV(a;)e aa: da;J H S(x)AT(a;)J = 2 ( y ^(^A^^e^dxJ Aq /2 , 

B(x)N(x)e Xox dx < 4A < 4B M . (24) 



we obtain that 

fOO 



/ 

JO 



'0 

Then, we can iterate on a and choose any a < B m + Ao < 2Ao while keeping the right hand side of 
(|23p bounded. This proves the L 1 statement of (p2l . 

It remains to use (|17[) again to deduce, thanks to the chain rule, that -^[e ax N(x)] € L 1 (M + ), and 
the result follows. □ 

With these results and methods of proof, we can now state the main new result of this section. 
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Theorem 3.2 Under the assumption U6\) , the map B \— > (Ao,iV), from L°°(R + ) into [B m ,BM] x 
L 1 n L°°(R + ) is 

(i) Continuous in B under the weak-* topology of L°°(M + ), 

(ii) Locally Lipschitz continuous in B under the strong topology o/L 2 (M + ) into L 2 (M + ), namely \2b\ ) 
and ( E3P below hold, for m large enough, 

(Hi) Class C 1 in the spaces of statement (ii). 

In statement (i) the weak-* topology of L°°(R + ) is natural from the assumption (|16p . For estimates 
(ii) and (iii), we have chosen the L 2 space which is used later on for the inverse problem where it plays 
an essential role. 



Proof. First of all, we introduce some notations that we will use throughout the proof. For two 
functions B, B satisfying (|16p . we set 

5B = B-B, A = \\B - B\\ L 2 iR+) . (25) 

For the corresponding solutions we define the differences 

SN = N-N, 5X = X -X . 

(i) This first statement follows thanks to the strong uniform estimates on N in Theorem 13.11 Indeed, 
consider a sequence B n that weak-* converges to B and the corresponding solutions N n to (|17p . The- 
orem [3J] provides us with uniform bounds on N n and thus with uniform bounds on • This means 
that the sequence N n is relatively strongly compact. Therefore, it converges to the unique N solution 
of Equation (|17p for the limit B. 



(ii) Firstly, we estimate the difference between the eigenvalues. As a consequence of (I17p for N 
(multiplied by (p) and (fT8|) for <p (multiplied by N), we obtain 

5\Jn^ = J 5Bip(^)N + J 5BipN. 

Therefore, according to the a priori estimates of Theorem 13.11 on ip and on N, 

\8X\ <C(B,B)A, (26) 
with C(B,B) = \\^\\ LOO(R+) \\(l + x)N\\L^+)/jNcp. 

Secondly, we estimate SN. Writing the difference of the solutions to the cell-division equation, we 
have 

-9-6 N + (A + B(x)) 5N = AB{2x)5N(2x) + SR(x), (27) 
ox 

with 

SR(x) = 45B{2x) N(2x) + (5X + 5B(x))N(x). (28) 
We have, since N has exponential decay at infinity using again Theorem 13.21 and for all m € N, 

||^(1 + x m )|| L2(M+) <C(B,B) A, (29) 
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for another constant C still controlled by B m an ABm- Also, since tp has at most linear growth at 
infinity, we know that 

/■oo 

\ \SR(x)\(p(x)dx < C(B,B) A. (30) 
Jo 

It is also useful to recall that the Equation (|27[) is spectral (N is the Krein-Rutman eigenvector) 
and one readily checks the solvability and the uniqueness conditions 



<p5R = 0, / 5N = Q. (31) 
Jr+ 

Our local lipschitz regularity result then follows from the two next lemmas 

Lemma 3.3 (Entropy dissipation) Assume that Equation M6\) holds and that N and <p are the 

solutions of (T7\ ) and (TE\). Then, the solution 5N to (27\ ), with the conditions [31\) , satisfies for all 
convex function H : R —* K, 

coo r _ 

dx 



roc r 

= r H '(f$) 6R W dx. 



(32) 



This lemma is a variant of the generalized relative entropy inequalities introduced in [171 [13] . Its 
proof will be presented in the Appendix^] We point out, however, the main difficulty that arises here: 
Since and N do not have necessarily the same exponential decay at infinity, one cannot use the 
natural weighted quantities J* °° ipN\j^\ p which are not necessarily finite. Therefore, we are sometimes 
forced to work in L . This type of entropy inequality is related to Poincare type inequalities which 
would be a way to conclude the proof at this stage. This area is still a very active field of research for 
so-called hypocoercive operators. See for example |18| . Here, it does not seem to follow from standard 
methods and we prefer to prove directly the following consequence of Lemma (|3.3p 

Lemma 3.4 (Spectral gap) There is a u > such that 

u\\8N\\ L 2 {R+) < \\5R (1 + x m )\\ L2(R+) < CA, (33) 
where the parameter m (large enough) is chosen so that Xq > ■ 
We again refer to the Appendix for a proof. 
This Lemma expresses the Lipschitz regularity stated in point (ii). 

(iii) The third claim is now a general fact on the Frechet derivatives of mapping with quadratic non- 
linearities. Whenever it is Lipschitz continuous on a subset, it is C 1 by usual algebraic manipulations. 
This concludes the proof of Theorem 13.21 □ 



8 



4 The Inverse Problem and its Regularization 



The inverse problem consists in finding B from the knowledge of the population growth rate and the 
cell density (Ao, N). If one could guarantee that the measurement N is very smooth one could directly 
consider to solve the equation on B, equivalent to the cell-division Equation (|17p written with y = 2x, 



4B(y)N(y) = B( y -)N( y -) + X N( y -)+2^- 



dy v 2 ; 



y >o. 



(34) 



This is a well-posed equation on B as long as iV satisfies regularity properties such as J^-^(f) ^ 
for some p > 1. But our interest lies on the fact that fluctuations on TV" make it a mere LP function 
because, given a set of measurements on N, one does not have a way of controlling the precision of 
the measurements on ^-N . This motivates working with a regularized problem 



' af y {B a N) + AB a {y)N{y) = B a (|)iV(|) + A iV(|) + 2^iV(|), y > 0, 
, B a N(0) = 0, 



(35) 



where < a < 1 is a small parameter adapted to the level of noise as explained in Section [2] for the 
toy model. Notice that this is still a well-posed problem in appropriate function spaces when written 
with a general source F(y) 



■ a-^(B Q N) + AB a {y)N(y) = B a {%)N($) + F(y), y > 0, 
, B a N(0) = 0, 



(36) 



Before we analyze the issue of estimating the inverse problem, let us explain why this is indeed a 
well-posed problem and for simplicity we restrict ourselves to L 2 . In particular, we did not impose 
boundary conditions because we expect that the point x = is characteristic since the cell density 
satisfies 



N G H 2 (R + ), N(x) > for x > 0, N(0) = ^N(0) = 0. 



(37) 



Theorem 4.1 (Strong stability) Assume (37) and that F £ L 2 . Then, there is a unique solution 
to ME) such that BN € H 1 with 



and 



a\B a N(y)\ 2 + J \B a N\ 2 <C j \F\ 2 , Vy > 0, 
a 2 J \^-B a N\ 2 <C f \F\ 2 , 

f\±B a N\><cf\±F\> ifF(0) = 0. 
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Proof. We drop the index a in order to simplify the notation. This model is a variant of the cell 
division equation and the Cauchy-Lipschitz theory applies and gives the existence for small y. The 
important point is to establish the a priori bounds. For the first set of estimates, we multiply the 
Equation (|36p by BN and integrate from to y. This yields 



%\BN{y)\ 2 + 4jv\BN\ 2 (s)ds 

= f V (BN)(z) (BN)£)dz + [ V (BN)(z)F(z)dz 
Jo 1 Jo 



<I f V \BN\ 2 + - l V \BN(-)\ 2 dz + - [ V \BN\ 2 + - r \F(z)\ 2 dz. 
~ 2 Jo ' ' 2 Jo 1 V 2j 1 1 2j 1 { M 

Altogether, this gives the two bounds claimed in the first statement above. 

The second estimate, on a J" \ J^BN\ 2 , is then a consequence of the first two using Equation ([36 
once again. 

For the third inequality, we differentiate the equation and write Q = J^BN. We find 
Therefore, multiplying by Q and integrating by parts gives 

/* OO 1 f OO f OO C\ 

l] Q 2 <^Q\0) + -J q Q{y)Q{\)dy + J Q(y) — F(y)dy, 



and thus 



Altogether this gives 

11 

T 



fV<?Q 2 (o) + Tl^l 2 - 



It remains to notice that from the Equation ([36]) on i?iV and the condition N(0) = we have also 
aQ(0) = a-§^BN(0) = F(0) and the result is proved. □ 

The first question we want to answer is how much this regularization on B differs from its exact 
value 

Theorem 4.2 (Consistency) Assume (37\ ), then the solutions to [34\ ) and (35\) satisfy 

\B„ - B\ 2 N < Co? f \\^-N\ 2 + \^N\ 2 

ay oy L 
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Proof. We write 5B = B a — B and we have 



a— (6B N) + A5B{y) N(y) = <5£(|) iV(|) - a^-(BN). 

This is again an equation of the form (|36|) with source term 

F=-a^-(BN) 
dy 

and thus we can apply the first estimate of Theorem 14.11 and obtain 

2 



J \SB N\ 2 < Co? j 



Next we apply the second estimate of Theorem 14.11 to (f34"|) with F = XqN + j^N and obtain 



Altogether these two estimates give the announced inequality. □ 



5 Recovering Estimates for the Regularized Inverse Problem 

Consider now a smooth solution (Xq,N) to Equation (|18|) corresponding to division rate B, that we 
want to recover from a noisy measurement iV £ with 



\N -N E \ 2 < e 2 , 



N_<N £ < iV 4 



(38) 



where N± are two smooth functions that can serve as a filter for unrealistic data. We assume their 
behavior contains possible properties of the true solution and serves to assert the boundary condition 
is well defined with -/V + (0) = iV + (oo) = 0. Notice N_ = is a possibility. The regularized method is 
aimed at furnishing an approximation B e ^ a of the exact coefficient B through Equation ()35p with N £ 
in place of N, 



a^(B £:a N £ ) + 4B £!a (y)N £ (y) = S £iQ (|)iV e (|) + A iV(f) + 2^iV £ (|), y > 0, 



(39) 



(B e>a N e )(0) = 0. 

We prove the following error estimate on the coefficient recovered by the above regularization procedure 



Theorem 5.1 (Convergence rate) Assume (37\ ) and [3B\) . then the solutions B to (&4\ ) and B £ a 
to \39\) . satisfy the error estimate, 



\Be,a - B\\ L 2( N 2 dx -) < Ca\\N\\ H 2 + 



C+\\B\ 



a 



N £ - N 



L 2 ■ 



11 



This theorem relies on a first estimate which expresses weak stability and improves Theorem 14,11 
using the special structure in the RHS of Equation (|39p . Namely, we have 



Proposition 5.2 (Weak stability) Assume |37[ ) and Hflfy . given two solutions B a to (fff5J) and IW\ 

satisfy 



\B £ , a N £ - B a N\ 2 < ^\\N £ - N\\\ 2 . 



Proof. We set Q = B £ ^N £ ^ a - BN, R = N £ - N which satisfies 

a-^Q + 4Q(y) = Q(|) + X R + 2^fl(|) . 

Then, we obtain after multiplication by Q and integration 

^[Q 2 (y)-Q 2 (0)]+4 [ V \Q(z)\ 2 dz = r Q ( z )Qd)dz + 
1 Jo Jo 1 

A J QR-2 J R(i)—Q + 2R(^)Q(y)-2QR(0). 

Using again the equation on Q (and the Cauchy-Schwarz inequality) we find 

f IQ 2 (V) - Q 2 (0)] +2 J* \Q(z)\ 2 dz < Cj Q y \R\ 2 + 2Q{y)R[l) - 2QR(0) 

%R(%)-4Q(z)]dy. 

Therefore 

^\Q(y)--R( V -)\ 2 + f y \Q(z)fdz<^\Q(0)-*R(0)f + C ( V \R\\ 
2 a 2 J 2 a a 1 J 

Notice that even though R is not defined point-wise, Q{y) - is well defined (as a H 1 function) 

and it vanishes due to (f38l) , □ 



We are now ready to prove our main result, Theorem 15.11 We write 
\\B £ , a - B\\ 2 L2{N!dx) < 2 / \B £:a N £ - BN\ 2 + 2 / \BN £ - BN\ 2 

< 4/ \B EiCt N E - B a N\ 2 + 4/ \ B a N - BN\ 2 + 2 / \BN £ - BN\ 2 . 

The first term is controlled thanks to Proposition 15.21 and gives the second error term in the estimate 
of Theorem 15. 11 The second term is controlled because of Theorem 14.21 and gives the first error term. 
The third term gives the last contribution once ||B||£oo has been factored out and the Theorem 15.11 is 
proved. □ 
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6 Conclusion and Suggestions for Further Research 



We have proposed a non-standard regularization method for the Inverse Problem associated with the 
cell-division equations. It is based on a first order operator that leads us to solving a new equation 
whose structure is very close to the equation itself (but for a different unknown). This allowed us to 
develop a consistency and convergence analysis that can be used in practice because the outcome are 
standard inequalities taking into account the noise in the measured data. 

Our theoretical approach has been focused on the equal mitosis, when mother cells divide into two 
equal daughter cells. A first extension of the method would be to deal with more general division 
equations as 



t n(t,x = 0) = 0,t > 0. 

The function (5 represents now the repartition of daughter cell sizes and a natural inverse problem 
would be to recover the division rate B from measurements of the cell density N once /3 is known. In 
principle, our method can extended but technical estimates have to be reformulated in this new context. 

The numerical validation of the procedure proposed here is straightforward because it relies on 
a regularized equation that has a standard form and numerical methods are available. As far as 
validation on real data is concerned, let us point out that present experimental devices allow us to 
measure not only the size repartitions in a cellular culture, but also the molecular content of certain 
representative proteins. Mathematical models describing the division process based on molecular 
contents are being developped. To our knowledge, identification of coefficients in this context is still 
largely open. 

Acknowledgments 

J.P.Z. acknowledges financial support from CNPq, grants 302161/2003-1 and 474085/2003-1. This 
research was made possible through the Brazil-France cooperation program, which is gratefully ac- 
knowledged by both authors. 

A Appendix 

This appendix is devoted to the proof of two technical lemmas that were used in the Section [3] to 
obtain Lipschitz regularity of the solution to the eigenproblem for the cell division equation. 

A.l Generalized Relative Entropy. Proof of Lemma 13. 3L 

Using ([2?]) and (fT7|) . we have 




(t,x) + £[g(x)n(t,x)]+B(x)n(t,x) = 2/ °° B(y)^)n(t,y) 



x>0,t>0 



d SN(x) 
dx N(x) 



AB(2x) 



N{2x) ,5N(2x) 5N(x) 
N(x) L N(2x) ~ N(x) 



5R{x) 
~N(x)' 



(40) 
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and thus, for H(-) a convex function, 

d nf 6N(x) , N(2x) ,, 6N(x) v 6N(2x) 5N(x) n , r SN(x)^ 5R(x) 



dx v N(x) > v ' N(x) v N(x) ' 1 N(2x) N(x) J v N(x) 1 N(x 
On the other hand, combining Equations (|17|) and (|18|) . we have 

— (N(x)<p(x)) = 4ip(x)B(2x)N(2x) - 2N(x)B(x)<p£). 

(JJL — ■ 

Therefore, 

I [N(xMx)H(^)] = AB(2x)N(2xMx)H>(™$) [fg - ™§>] 

+ [A^(x)B(2x)N(2x) - 2N{x)B{x)y{l)]H{*0±) 
+H>(§$) SR(x) V (x). 

After integration in x we arrive at 

f°° . r sr., \ i\t / .„,,6N(x), r 6N(2x) SN(x) 1 , 
M ,)B(2,)i V (2x)/ f '(^)[^i-^i ; i]< i x 



+ fM I )B( 2l )A f ,2,)[«(^|)- H (^g))] d , 



+ / H'( S -^±)6R(x)<p(x)dx. 



o 



Once the terms are reorganized, this is exactly the statement of Lemma 13.31 □ 



A. 2 A Poincare-Type Inequality. Proof of Lemma 13.41 

We argue by contradiction and assume that for a sequence B k — > B, the corresponding v k vanishes. 
Then, there is a family of 5R k , 5N k that satisfy 

v k \\6N k \\ L2(M+) > \\6R k (l + x m )|| L2(M+) . (41) 

After re-normalizing (by multiplication), we can always assume that 

l|<5iV fe || L 2 {R+) = 1, \\SR k {l + x m )|| i2(R+) 0. 

Then, our proof consists in several steps: First, we prove compactness and then pass to the limit, next 
we use the generalized relative entropy to identify the limiting solution and finally, prove a contradic- 
tion. 

Compactness. Notice that the Equation (|27p automatically implies that (recall (p > is smooth and 
sublinear at infinity), that SN k is bounded in H 1 (W + ). This provides us with local compactness. It 
remains to obtain a control at infinity; we also notice that 

\^ N k? + (Ao + B(x)) (5N k ) 2 = 4B(2x)5N k (2x) 5N k (x) + 5R k (x) 5N k (x) , 
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which, after integration in x with the weight x m gives for all a > 

-f Jr+ z™- 1 ^) 2 + / R+ x m (A + (5iV fe ) 2 < 4/ R+ x m B(2x)[a5N k (2x) 2 + ±5N k (x) 2 } 



+^/ R+ x-(5iV fc )2/ R+ x ™(6R k ) 2 . 

Therefore, we also obtain 

J R+ x m (X + B(x)) (5N k ) 2 < S J R+ x — 1(5^)2 + a 2~™+ l f R+ x m B{x)5N k {x) 2 



With a = 2 m_1 , we conclude that 

A / R+ *"W fc ) 2 < § J R+ x" 1 " 1 ^) 2 + / R+ ^ m (<5iV fc ) 2 + ^/ R+ x™(5Ay 2 / R+ * m («4) 2 

< C/ R+ (<WV fc ) 2 + (4f + ) / R+ x m (<W fc ) 2 + ^/f R+ x™(5N k yf R+ x™(5R k ) 2 . 
It remains to choose m large enough so that Aq > 2^M. anc [ we conclude on a uniform control 



x m (8N k ) 2 < C 



1+ / x m (5R k f 



(42) 



that concludes the global compactness in L 2 (M + ) of the sequence SN k . 



Limit. Therefore, due to the above facts and using Lemma 13, 3\ we may pass to the strong limit and 
obtain 

5N k -» n G x m n G L 2 (M+), V?n > 0, 



<5n = 0. 



5 r 
—n + (X + B(x))n = AB(2x)n(2x), 
ox J R 

Notice that, thanks to the moment estimate, we also have for all m > 0, 

(1 + x m )n G and / \n\tp < oo. 

This allows us to give a meaning to the vanishing integral. This follows from the Cauchy-Schwarz 
inequality 

ft r dr \ 1/2 

(l+x m )\n\<( (l+x m ) 2 (l + \x 2 \)\n\ 2 (x)dx - = 



Application of the GRE. We can now choose the family of convex functions H (u) = (u — £)+ in the 
Generalized Relative Entropy of Lemma 13.31 with £ > a parameter. We find, because we now deal 
with 5R = 0, after reorganizing the terms in (I32|) . that 
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In fact, one needs to justify it by a preliminary truncation of the integral at x = and x ~ oo, we leave 
to the reader the corresponding analysis. Again, let us point out that in L 1 this is possible because 
the corresponding ^-derivative term 

is integrable (but for the square entropy we arrive at [ip(x)N(x) ( ) 2 ] which is not integrable!). 
Conclusion. From the previous step, we conclude that, for all x > 0, £ > 0, we have 

Sgn+( iv>)-^ =sgn+ (iv(^)- e )' 

which implies = jfi^j an d thus, tells us that the limit satisfies 

d n(x) 



dx N(x) 
'N and tl 

This gives the sought contradiction and concludes the proof of Lemma 13.41 □ 



From jffil = ^2x) ' we concm de that n = CN and the vanishing total integral of n implies that n = 0. 
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